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Objectives 


Accurately model the raw GNSS signal received during rendezvous and close proximity 
operations around large space structures 


Validate this model using experimental data 


Use this information to improve relative navigation accuracy and integrity 
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Discovery docking with ISS May 29, 1999 [1] 
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Discovery docking with ISS May 29, 1999 [1] 
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Background 


During docking, GNSS navigation of the chaser spacecraft is corrupted by GNSS signals 


reflected off the passive target spacecraft 


These reflected signals, i.e., multipath, could potentially be used as a source of information 


about the distance of the target from the receiver 
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Background (continued) 


A model has been presented previously [2] that uses electromagnetic (EM) ray tracing and a 
MATLAB signal simulator to generate raw, multipath-corrupted GNSS data from environment 


features 


Hubble Servicing Mission 4 (HSM4) is used as a case study, and space data collected during 


this mission is used here to validate the model 
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Background — Signal Model 


. Received direct signal (GPS C/A signal on the L1 carrier) 
y(t) = ad(t — Teode(t) )e(t — Teode(t)) cos(2m frit + O(t)) 
with carrier phase 
O(t) = —27 fpr Tearr(t) + 90 


- Reflected signals are received as delayed, attenuated replicas of the line of sight signal. For M 
reflected signals: 
M 


ai S- Amd(t — Teode,m(t))c(t — Tcodeym(t)) cos(27 frit + Om(t)) 


m=0 
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Background — Signal Model 


Received direct signal (GPS C/A signal on the L1 carrier) 
y(t) = ad(t — Teode(t) )e(t — Teode(t)) cos(2m frit + O(t)) 
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reflected signals: 
M 


ai S- Amd(t — Teode,m(t))c(t — Tcode,ym(t)) cos(27 frit + Om(t)) 


m=0 


- Characterize reflected signals in terms of relative amplitude, code delay, and carrier phase with 


respect to the line of sight signal: 
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Background — Signal Model 


. Received direct signal (GPS C/A signal on the L1 carrier) 
y(t) = ad(t — Teode(t) )e(t — Teode(t)) cos(2m frit + O(t)) 
with carrier phase 
O(t) = —27f rr Tearr(t) + 90 


- Reflected signals are received as delayed, attenuated replicas of the line of sight signal. For M 


reflected signals: 
M 


1 S- Amd(t — Teode,m(t))c(t — Tcodeym(t)) cos(27 frit + Om(t)) 


m=0 


- Characterize reflected signals in terms of relative amplitude, code delay, and carrier phase with 


respect to the line of sight signal: 


| i bm(t) = C(Teode,m(t) — Teode,o(t)) 
Am = | i a 
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Background — Signal Model 


. Received direct signal (GPS C/A signal on the L1 carrier) 
y(t) = ad(t — Teode(t) )e(t — Teode(t)) cos(2m frit + O(t)) 
with carrier phase 
O(t) = —2 fpr Tearr(t) + Oo 


- Reflected signals are received as delayed, attenuated replicas of the line of sight signal. For M 


reflected signals: 
M 


ai S- Amd(t — Teode,m(t))c(t — Tcodeym(t)) cos(27 frit + Om(t)) 


m=0 


- Characterize reflected signals in terms of relative amplitude, code delay, and carrier phase with 


respect to the line of sight signal: 
sm Boy OH = [8 din (t) = C(Teodesm(t) — Teateo(t)) 
mT Gro Vag Ym(t) = Go(t) -Om(t) 
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Background — Multipath Fading 


As the excess path length traveled by a reflected signal changes the relative carrier phase will 
cycle through phases that add to the line of sight signal constructively and destructively 


This effect is known as multipath fading, resulting in a characteristic oscillation of received power 
Multipath fading is driven by the carrier phase rate of change: 


Ym(t) = hen (t) 


This can be related to the rate of change of the relative code delay: 


dint) = 5— (2M) dnt 


carr 


Ignoring dispersion effects: 
1 
Um(t) © ao oe b(t) 


This relationship is used to validate the signal simulation - the simulated reflected signal path 


length rate of change is compared to the measured fading frequency. 
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Hubble Servicing Mission 4 


STS-125 was performed by the crew of the 
space shuttle Atlantis May 11-24, 2009; 
docked with Hubble on May 13" 


Relative Navigation Sensor (RNS) experiment 


¢ Tested vision processing algorithms that 
estimated Hubble's relative position and 
attitude using imagery from three 


cameras mounted in the shuttle bay 


¢ Two oppositely-polarized GPS antennas 
were mounted in the shuttle cargo bay 
(RHCP and LHCP) collecting hours of 
sampled, intermediate frequency GPS 


data 
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Atlantis docking with Hubble May 13, 2009 [3] 
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HSM4 - Antennas 


- Sensor Systems RCP and LCP antennas flown to 
receive direct and reflected signals respectively 


- Separability of direct and reflected signals depends 
on polarization reversal at reflection and antenna 


cross-pole discrimination 


- Antennas were measured in the Goddard 
ElectroMagnetic Anechoic Chamber (GEMAC) 


¢ Attenuation of LCP signal at RCP antenna 
boresight is 12 dB 


¢ Attenuation of RCP signal at LCP antenna 
boresight is 7 dB 


¢ Azimuthal variation due to asymmetric 


mounting plate 
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Antenna testing in GEMAC with antenna mount-plate reconstruction 
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HSM4 — Antennas (continued } 


RCP Gain 
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Simulation 


Simulation presented previously [2]: 


¢ Dynamics of Hubble and GPS 
satellites reproduced in the 
WinProp© electromagnetic ray 
tracing software using a CAD 
model of Hubble 


¢ Modified GSFC Siggen was used 
to produce GPS data from the ray- 
tracing results and a custom 
MATLAB software receiver was 
used to track simulated and HSM4 


data 


National Aeronautics and Space Administration 


16:31:00 


16:31:10 


16:31:20 


16:31:30 


16:31:40 


16:31:50 


PRN 28 interactions 


16:28:00 


16:28:20 ae ae 
| 

16:28:30 a ar 
| | | 


16:28:50 


| | 


\ 


PRN 17 interactions 


ION GNSS+ 2016 Portland, OR 


simulation (continued) 


Ray tracing output consists of relative path length and electric field strength for each ray 
¢ Simultaneous arrival of nearly identical rays resolved through curve fitting 


* 165 second segment considered here: May 13", 2009, 16:25:27-16:28:12 UTC (inter- 
vehicle range ~85 meters) 
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Results 


Simulated parameters of reflected signals for four PRNs (FFT resolution = 0.0061 Hz) 
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PRN | do [m] | 6 ea) a2, /a2 | wh [Hz] 
1 | 70.80 0.2595 
26 | 1675 | 044 | 0.66 | 0.0230 
27 | 1632 0.0246 
30 73.84 5.73 0.38 0.3012 
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Results — PRN 1 (simulation) 


Simulated parameters of reflected signals for four PRNs (FFT resolution = 0.0061 Hz) 


9 x 40° PRNO1 Prompt Correlator Power 


RCP corr. power 
& ol Oo 


w 


16:26 


National Aeronautics and Space Administration 


16:27 
time [UTC] 


7080 | -494 | 045 


16:28 


0.2595 


PRNO1 Prompt Correlator PSD 


eal 
a 
oO 


0.2595 Hz 


ak aks 
ie) Py 
oO oO 


» 


Bare 
Nh 
oO 


RCP correlation PSD [dB] 
ry 


0.5 1 15 2 2a 
frequency [Hz] 


ION GNSS+ 2016 Portland, OR 


Results — PRN 1 (experimental) 


- Simulated parameters of reflected signals for four PRNs (FFT resolution = 0.0061 Hz) 


RCP corr. power 


LCP corr. power 


16:26 16:27 
time [UTC] 


0 x 10° PRNO1 Prompt Correlator Power 


“| 


2 


hd hk i eer y | 
16:26 16:27 
time [UTC] 


National Aeronautics and Space Administration 


16:28 


16:28 


PRNO1 Prompt Correlator PSD 


130 

a —RCP power 
% 120 — 20 ms noncoherent avg 
Aa 

ep) 

oO 

c 

2 

% 

© 

8 

o 

oO 

fe 


0 0.5 1 1.5 2 2.5 3 


frequency [Hz] 

430 PRNO1 Prompt Correlator PSD 
— ——LCP power 
g 120 —— 20 ms noncoherent avg 
D 110 
= 0.2594 Hz 
o 
ro 90 
8 80 
Oo 
O 70 


frequency [Hz] 


ION GNSS+ 2016 Portland, OR 


23 


Results — PRN 30 (simulation) 


- Simulated parameters of reflected signals for four PRNs (FFT resolution = 0.0061 Hz) 
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Results — PRN 30 (experimental) 


- Simulated parameters of reflected signals for four PRNs (FFT resolution = 0.0061 Hz) 
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Results — PRN 26 (simulation) 


Simulated parameters of reflected signals for four PRNs (FFT resolution = 0.0061 Hz) 
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Results — PRN 26 (experimental) 


- Simulated parameters of reflected signals for four PRNs (FFT resolution = 0.0061 Hz) 
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Conclusions 


Multipath fading confirms the presence of reflected signals in the HSM4 data 


General agreement of these features with the effects of simulated Hubble-reflected signals 
provides validation that the multipath model correctly represents the GNSS reflections received 


during docking 
Ray tracing results support the hypothesis that there are reflections off Hubble of sufficient 
strength and duration to be useful for relative navigation 
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Future Work 


- Resolve azimuth uncertainty and incorporate complete antenna gain patterns 
- Extract range information from reflected signals, as explored in other work [4] 


¢ Estimation of multipath parameters from code correlation deformation — relies on a wide 
receiver bandwidth (HSM4 data bandwidth = 2 MHz) 


¢ Independent tracking of reflected signals — relies on reflected signal isolation (e.g., LCP 


cross-pole discrimination > 10 dB) 


Apollo and Soyuz docking 1975 [5] 
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Backup — Signal Model 


If we let O(t) = —2rfriToarr(t) + Oo, then 
y(t) = ad(t — Teode(t))e(t — Tonde(t)) cos(2x frit + O(t)). (5) 


Note that if Doppler is approximated as constant and fp = f;,:7(0)t/c is substituted, #(t) = —2x fx,r(0)/ce+ 
fot. The received signal can be rewritten by grouping f,; and fp, with the phase 6(t) = —2 f,,r(0)/e+ 


Oo: 
; y(t) = ad(t — teose(t))e(t — teoue(t)) 208(2a( fur + fo)t + 6(0)), © 
Teode(t) = —— — (r(0) +#(0)t+ 0) 
Tearr(t) = — (r(0) + F(0)t + O) 
Yin(t) = Oo(t) — Om(t) = —2% fx (Fearr,o(t) — Toarrm(t)). 


edie 


¥in(t) = — —2r fir 


(xO - sat _, Pan( 0). 


Vearr 


vat) > in(t) = 


bm(t) = CC Feade yn (t) — Feode,o(t)) = = (Fm(t) — Fo(t)) = zfamlO). 
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Backup — PRN 27 (simulation) 


Simulated parameters of reflected signals for four PRNs (FFT resolution = 0.0061 Hz) 


PRN | 60 [m] | 6 [cm/s] | a2,/a2 | w [Hz] 
27 | 16.32 0.47 0.67 0.0246 


~4 periods / 165 sec = 0.0242 
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Backup — PRN 27 (experimental) 


Simulated parameters of reflected signals for four PRNs (FFT resolution = 0.0061 Hz) 


PRN | 60 [m] | 6 [cm/s] | a2,/a2 | w [Hz] 
27 16.32 0.47 0.67 0.0246 


Note: LCP tracking of HSM4 data failed for PRN 27 and is not shown 
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